Can cosmic strangelets reach the earth ? 
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The mechanism for the propagation of strangelets with low baryon number through the atmo- 
sphere of the Earth has been explored. It has been shown that, under suitable initial conditions, 
such strangelets may indeed reach depths near mountain altitudes with mass numbers and charges 
close to the observed values in cosmic ray experiments. 
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The existence of Strange Quark Matter (SQM), containing a large amount of strangeness had been postulated by 
I various authors quite a few years ago. In a seminal work in 1984, Witten ^ proposed that SQM with roughly equal 

■ numbers of up, down and strange quarks could be the true ground state of Quantum Chromodynamics (QCD), the 
' accepted theory of strong interactions. While only SQM with very large baryon numbers were initially thought to be 

favorable (in terms of stability), later calculations have shown [^-^ that small lumps of SQM can also be stable. The 
^ I occurrence of stable (or metastable) lumps of SQM, referred to in the literature as strangelets, would lead to many 
, ^ I rich consequences; for a recent review, see [0. 

I Strangelets may arise from various scenarios; they could be formed in highly energetic nuclear collisions associated 
with the formation of quark-gluon plasma |^J^, or they might be of cosmological origin, as remnants of the cosmic 
QCD phase transition Collisions of strange stars could also lead to the formation of strangelets which could 
I contribute to the cosmic ray flux 0. In heavy ion collisions, it is thought that strangelets with atomic number A 
J> [ upto 20 - 30 may be formed , the stability of which depend rather sensitively on the parameter values (like the Bag 
' constant) and an underlying shell- like structure. For larger strangelets ( ^4 > 40), the stability appears to be more 
robust I^J^. We confine our attention in this work to such larger strangelets, which may not be readily formed in 
r heavy ion collisions in the laboratory but could be of cosmic origin. A discerning property of such strangelets would 
[ be an unusual charge to mass ratio ( Z/A<^ 1) 0. 

. The obvious place to look for such strangelets would be in the cosmic ray flux. In this context, it may be recalled 

■ that there have been intermittent reports in the literature p|-p^ about the detection of exotic cosmic ray events, with 
(— I unusually low charge to mass ratios; some of these events are tabulated in Table 1. Although it appears natural to 

identify these events with strangelets, no consensus has yet emerged, primarily because of the ambiguities associated 
^ I with the mechanism of propagation of strangelets through the terrestrial atmosphere. For example, if a strangelet 
^ ' arriving at the top of the atmosphere has a baryon number A ~ 1000, there would be a serious problem with its 
. penetrability through the atmosphere, as the exotic events are observed at quite low altitudes. One could assume that 

■ their geometric cross sections are very small. Alternately, one could conjecture, a la Wilk et al |l^-[l6|| and others [pT| , 
that although the initial mass of the strangelet is very large, it decreases rapidly due to collisions with air molecules, 
until the mass reaches a critical value merit below which the strangelet simply evaporates into neutrons. 

The difficulties associated with this kind of interpretation are twofold. Firstly, one has to take account of the fact 
that, unlike ordinary nuclear fragments which tend to break up in collisions, strangelets can become more strongly 
bound if they absorb matter [Q . Secondly, since a strangelet has a net electric charge, it experiences an ever increasing 
geomagnetic field, which considerably lengthens its path before reaching a certain altitude. This implies many more 
interactions with the nuclei of the atmospheric atoms, as a result of which the strangelet would "evaporate" much 
before the desired depth is reached. 

These difficulties can be naturally overcome in a different scenario, proposed recently by the present authors [ p^ , 
in which the stability of the strangelet plays a very important role. In this model, an initially small strangelet, during 
its travel through the Earth's atmosphere, picks up mass, rather than lose it, from the atmospheric atoms. Such a 
situation may prevail unless the propagation velocity of the strangelet through the atmosphere is so high that in a 
collision with the atmospheric nucleons, the excitation energy would exceed the binding energy. We have estimated 
that for our case, where the initial A is larger than 40, this upper limit on the velocity comes out to be above 0.7c. 
(We disregard the possibility of fission-like fragmentation of the strangelets.) The equation governing the rate of 
change of mass with respect to distance traveled is given by, 

drus _ f X nin . . 

dh ^ X ^ ' 

and the equation of motion reads 
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In the above equation, tUs and v represent the instantaneous mass and the velocity of the strangelet, q is the charge 
and A is the mean free path of the strangelet in the atmosphere. The factor / determines the fraction of neutrons that 
are actually absorbed out of the incident neutrons (to„ ). In this case, A is both a function of h (which determines 
the density of air molecules) and (which is related to the interaction cross section). The initial velocity has to be 
bigger than a threshold value, so that a strangelet of a given initial mass and charge can arrive at an altitude 25 km 
from the sea level , surmounting the geomagnetic barrier. The upper limit of 25 km is chosen primarily to economise 
on the computation time and is a fortiori justified since the density of the atmosphere above this height is almost 
negligible for our purpose. The variation of atmospheric density with height has been described by a parametric fit 
to the data given in the standard reference of Kaye and Laby . 

According to the above analysis, a strangelet with an initial mass of ~ 64 amu and charge ~ 2 evolves to a 
mass ^ 340 amu or so, by the end of its journey, an altitude ~ 3.6 km above the sea level (typically the height 
of a north east Himalayan peak in India, like Sandakphu , at a geomagnetic latitude ~ 30^ N). This mass is quite 
close to the few available data (see Table 1) and seems to support the interpretation that exotic cosmic ray events 
with very small Z/A ratios could result from SQM droplets. However, it was assumed in that only neutrons 
are absorbed preferentially over the protons from the nuclei of the atmospheric atoms (i.e. charge of the strangelet 
remains constant), the protons being coulomb repelled. It should nonetheless be realized that in the earlier phase of 
the journey, when the relative velocity between the strangelet and the air molecule is large, some protons will indeed 
be absorbed, albeit with a lower cross section than that for neutron capture. As the strangelet builds up in mass as 
well as in charge, the coulomb barrier at the surface of the strangelet gets steeper and the relative velocity also gets 
further reduced. This will slow down the charge transfer process and ultimately inhibit it. Also, one cannot avoid the 
issue of loss of energy of the strangelet through ionisation of the surrounding media. As we shall see, the ionisation 
losses, which become quite significant at comparatively low altitudes, actually provides a lower limit to the height at 
which the strangelets can be detected successfully. 

In this letter, we therefore try to explore the consequences of the absorption of protons by the strangelets in 
course of their journey through the terrestrial atmosphere in a relativistic setting. The equation of motion (H) can be 
generalized to a relativistic form in a straightforward manner: 

7m,- = -m,g + q{v x B) - + ^)- ^r^sV- - -j=-v (3) 

where 7 is the Lorentz factor. The third term takes care of the deceleration of the strangelet due to the absorption 
of neutron as well as protons, where the proton absorption term is related to the neutron absorption term as 

dms^ _ Op_ dm^^ ^ ^ dm^^ 
dt cr„ dt ^" dt 

where ap and (t„ are the cross sections for neutron and proton absorption, respectively. Treating, classically, the 
proton of energy E as a, free charged particle of unit charge in the repulsive coulomb field of the strangelet, we can 
easily estimate the minimum separation r,„i„ along the trajectory to be given by 

1^ + C/(.„.„) = E 

where U{r) represents the potential energy of the proton due to the coulomb field of the strangelet; Vo is the relative 
speed with which the N2 nuclei (and hence, its constituent protons) approach the strangelet and h is the impact 
parameter. Assuming that charge transfer can take place when rmin < Rs (the radius of the strangelet), the corre- 
sponding value of b{= he) is &^ = i?s^(l — U{Rs)/E), so that the proton capture cross section ( ct^ ) by a strangelet 



of atomic number Zs is ap = nb^ = ttR^ 
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In contrast, the scattering cross section for neutrons (cr„) is just 7r(r„ + i?s)^ and hence. 



{rr,+Rsf\ EAireoR. 

Finally, the last term of equation (||) accounts for the ionisation loss. The expression for f{v) is given by |^ 



2 



/ (^) ^ _^ = ^'"'"^"f ( (6) 

dx ineQlJleV^ \ bmm J 

Here, n represents the number density of the atmospheric atoms at a particular altitude, Zmed is the number of 
electrons per atom of N2 which can be ionised, me is the mass of the electron and bmax and 6mm are the maximum 
and minimum values of the impact parameter. At large velocities, expression (^) reduces to, with / denoting the 
average ionising energy. 

However, when the velocity of the strangelet falls below a critical value v < 2ZsVo{ vq — 2.2 x lO^m/s is the 
speed of the electron in the first Bohr orbit), electron capture becomes significant which can be accounted for by the 

replacement Zs Zl ^ io[|2|. 

Equation was solved by the 4*'* order Runge-Kutta method with different sets of initial mass, charge and p. It 
may be mentioned at this point that the first term in eqn is not important in magnitude, as is to be expected. 
We have nonetheless included it for numerical stability. This serves to define the downward vertical direction in the 
vector algorithm, especially for very small initial velocities. 

In figure 1, we have plotted final masses (for initial masses 42, 54, 60 and 64 amu ) with initial /3 for a fixed 
initial charge 2. This graph shows the following interesting feature; the final value of the mass decreases at first with 
increasing values of the initial /3 and then begins to increase again after a critical value of /3 is reached. This feature, 
although not apparent from the curve corresponding to the initial mass M = 64 amu, clearly reveals itself for lower 
initial masses. From the same figure, it can also be inferred that the value of /3 where the 'dip' occurs shifts to the 
left with increasing values of the initial mass. Although mathematically delicate ( it can be seen from eqn(^) that a 
higher value of speed leads to an increasing value of the mass increment, which in turn slows down the particle), a 
qualitative explanation of this feature might be given as follows. One can think of the total region through which the 
strangelet travels being divided into two distinct subregions. In subregion I, corresponding to higher altitudes, the 
number of atmospheric particles is small, while this number is considerably larger in subregion II, corresponding to 
lower altitudes. For small initial mass (smaller size), the strangelet has a greater chance to escape subregion I if /3 is 
higher, so that it will pick up lesser mass from this region. On the other hand, if /3 is very high, the volume that the 
strangelet sees will be contracted (the twisted tube through which it travels will be constricted) as a result of which 
it will interact with a greater number of atmospheric particles whence it will pick up a larger number of nucleons. It 
is clear that for an initially bigger (more massive) strangelet, this critical value of (3 will be lower, as it will be able 
to sweep through a larger number of atmospheric particles right from the start. 

Let us now consider a representative set of data with initial mass 64 amu and charge 2 for detailed discussion. The 
results for /3o = 0.6 are shown in figures 2 and 3, where the variation of speed (/3) and the energy of the strangelet 
with altitude are depicted. The sharp change seen at ~ 13 km corresponds to the onset of electron capture, which 
is handled phenomenologically through the effective Zg. The insets of figures 2 and 3 show a zoomed- up view of the 
respective quantities near the endpoint of the journey. It is apparent from the figures that the ionisation term reduces 
the overall energy and speed considerably from the nondissipative situation [^8|. However, the zoomed-up insets in 
figs. 2 and 3 show that the strangelets may have enough energy to be detectable at an altitude of 3.6 km from the sea 
level. For example, for the values of the initial quantities rris^ and /3o shown here, the strangelet is left with a kinetic 
energy -8.5 MeV (corresponding to ^ = 2.35 MeV/mg/cm^ in a Solid State Nuclear Track Detector (SSNTD) like 
CR-39), which, although small, is just above the threshold of detection {^)crit ~ 1 — 2M eV j mg j err? for (5 < 10~^ 
in CR-39 for the present configuration. Below this height, the possibility of their detection with passive detectors like 
SSNTD reduces to almost zero. 

Table 2 lists the final values of the quantities mass, charge, /3, and the energy of the strangelet at the end of the 
journey for different initial velocities. A comparison between tables 1 and 2 shows that the final masses and charges 
are very similar to the ones found in cosmic ray events. 

In conclusion, we have presented a model for the propagation of cosmic strangelets of none-too-large size through 
the terrestrial atmosphere and shown that when proper account of charge and mass transfer as well as ionisation loss 
is taken, they may indeed reach mountain altitudes, so that a ground based large detector experiment would have a 
good chance of detecting them. 
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of India, New Delhi. 
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Charge 


Counter experiment [^ 


A ~ 350-450 


14 


Exotic Track Eol 


A ~ 460 


20 


Price's Event [|ll| 


A > 1000 


46 


Balloon Experiments il^,|l3[ 


A ~ 370 
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TABLE I. Mass and charge obtained from cosmic ray experiments 
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418.3 
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456.3 
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TABLE II. The final values, denoted with suffix I, are tabulated along with initial /3 (/3o) 
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FIG. 1. Variation of final masses with initial /? for different initial masses 
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FIG. 2. Variation of final (3 with altitude (a) for constant charge and without ionisation loss and (b) including proton 

absorption as well as ionisation loss 
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FIG. 3. Variation of kinetic energy with altitude 
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